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THE PUZZLING CASE OF 56 PEGASI: 

A FAST ROTATOR SEEN NEARLY FACE-ON? 


By A. Frankowski and A. Jorissen 
Institut d’Astronomie et d’Astrophysique 
Universite Libre de Bruxelles 


A spectroscopic orbit has been recently found by R. Griffin for the long known barium 
star 56 Peg, which is also a strong X-ray source. This short-period, low mass-function orbit 
raises several questions regarding the history of the system. In the present paper, we show 
that it is not easy to find an evolutionary history for 56 Peg which is consistent with the 
current component’s masses, unless one assumes that the giant is a relatively fast rotator 
(a few times 10 km s^ 1 ). The hypothesis of fast rotation allows us to explain some other 
peculiarities of this object as well. 

[...] when you have excluded the impossible, 
whatever remains, however improbable, 
must be the truth. 

Sir Arthur Conan Doyle, 
The Adventure of the Beryl Coronet 


Introduction 

In paper 186 of R. Griffin’s series on Spectroscopic Binary Orbits from Photoelectric Radial Velocities 
appearing in this issue of the Magazine 1 , a 111.15 d orbit is found for 56 Peg, a rather luminous barium 
star which is as well a strong X-ray source. The star 56 Peg is also remarkable in having, despite its large 
radius (the object is classified as a bright giant, listed as K0.5 II by Keenan & McNeil 2 ; an interferometric 
measurement by van Belle et al. 12 yields a radius of 40±4.7 Rq for a uniform-disc approximation), the second 
shortest orbital period among barium stars, as revealed by Fig. Q This figure displays the eccentricity - 
period diagram of barium stars. Its mass function of 3.73 x 10~ 5 M 0 is the smallest among barium stars. The 
giant is of solar or slightly subsolar metallicity 3,4 : [Fe/H]= — 0.15±0.28. Older metallicity determinations 5,6 , 
in the range [Fe/H] = —0.21 to —0.01, are thus consistent with these recent values. For a thorough discussion 
of the object’s spectral type and abundances the reader is referred to the review by Griffin 1 . 

An absolute magnitude of My ~ —l m .3 ± 0.1 is quoted by Griffin for 56 Peg, based on literature values 
of the V magnitude and on its Hipparcos parallax' of 6.07 ± 0.67 milliarcsecond* (d = 165 ± 18 pc). 
Distance estimates based on the straight inversion of the parallax, as done here, are subject to statistical 
biases. In this specific case however, maximum likelihood estimates of the distance (based as well on the 
Hipparcos parallax, but providing unbiased distance estimates) yield distances consistent with the above 
value: Mennessier et al. 8 obtain d = 161 ± 16 pc while Famaey et al. 9 give d = 167 ± 18 pc. 

In what follows, wherever clarity requires to differentiate the giant component in 56 Peg from the system 
as a whole, we will use the name 56 Peg A for the giant. 

Location of 56 Peg A in the Hertzsprung-Russell diagram 

Several arguments presented in this note will be based on the location of 56 Peg A in the Hertzsprung- 
Russell (H-R) diagram. Therefore, we prefer to rely on consistent sets of magnitude and colour measurements 
from Tycho-2 10 and Hipparcos 7 . There are 193 Vt and Bt measurements available for 56 Peg, yielding 
< Vt > = 4.906 ± 0.022 and < Bt — Vt> — 1-53 ± 0.04, with the quoted uncertainty corresponding to the er 
of gaussians fitted to the Vt and ( B T — Vt) distributions (Hence, the standard error on the quoted mean is 
in fact smaller by the factor 193 1 / 2 = 13.9). The average uncertainty on single Vt and Bt measurements of 

* abbreviated ‘mas’ in the remainder of this paper 
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56 Peg is 0.021 and 0.038 mag, respectively, so that there is no indication from the Tycho-2 photometry that 
the star is variable (since the standard dispersions of both the Vp and (Bp — Vp) distributions simply reflect 
the measurement errors). Using Eqs. 1.3.29 and 1.3.33 of Volume 1 of the Hipparcos catalogue 7 to convert 
the Tycho-2 photometry to the Johnson system, one finally obtains V = 4.76 ±0.02 and B — V = 1.28±0.03. 
An analogous procedure using even more accurate measurements from Hipparcos Hp band leads to the same 
average value of V , but Hipparcos data do reveal small intrinsic photometric variability with a scatter of 
0.011 mag (Hipparcos catalogue field H46; see also the discussion in Griffin’s paper 1 ), which is below Tycho-2 
detection limit and is also sufficiently small not to pose a threat to our further reasoning. The obtained value 
for the V magnitude falls in the literature range V = 4.7 - 4.8 quoted by Griffin, and yields My = —1.32^'^g 
(where the quoted uncertainty combines the uncertainty on the parallax and on the photometry, but is in 
practice completely dominated by the parallax formal error). 

As the companion is much fainter that the giant and is not seen in the visible light (though it has been 
detected in the UV 1 ), we note that it is safe to ascribe the above photometry to the giant star. These values 
have been used in Fig. Q to locate 56 Peg A in the observational H-R diagram displaying stellar evolution 
tracks from Lejeune & Schaerer 11 , computed up to the first thermal pulse on the asymptotic giant branch 
(AGB). They reveal that in the solar-metallicity case (Z = 0.02), the giant star in the 56 Peg system has 
a mass of the order of 3-4 Mq, whereas the mass may be as low as 2 M 0 for the subsolar metallicity of 
Z = 0.008 (corresponding to [Fe/H] = —0.4). Admittedly, such estimates may depend upon the details of 
the complex theoretical modelling of late stages of stellar evolution. 

More importantly, as we will now show, 56 Peg A has not attained the thermally-pulsing asymptotic 
giant branch (TP-AGB) yet. Therefore, the nucleosynthesis associated with thermal pulses and the ensuing 
third dredge-ups cannot be invoked to account for the barium enrichment observed in the atmosphere of 
56 Peg A. 

From the observed parameters of 56 Peg A (B — V = 1.3 and My = —1.3) applying a bolometric 
correction one gets a luminosity of log L/Lq = 2.7. The luminosity at the onset of the TP-AGB has been 
estimated in the following way. In the solar-metallicity case (Z = 0.02; right panel of Fig.^J, the core-mass 
- luminosity relationship from Blocker 13 (as given by Eq. 1 of Herwig et al. 14 ) predicts logL = 3.5 at the 
start of the TP-AGB for a 3 M 0 star, when the core mass amounts to 0.54 M 0 according to Fig.0 Thus 
in the solar-metallicity case 56 Peg A is well below the luminosity defining the onset of the TP-AGB. The 
situation is similar for the Z = 0.008, 2.5 M 0 case (left panel of Fig. [2j : for a core mass of 0.52 M 0 at the 
onset of the TP-AGB (see Fig. 0. Blocker’s core-mass - luminosity relationship predits logL = 3.3, larger 
than the estimated stellar luminosity of logL = 2.7. The more complex core-mass - luminosity relationship 
of Wagenhuber & Groenewegen 15 yields basically the same results. Even considering TP-AGB stars in 
the luminosity dips following a thermal pulse (rather than in the quiescent H-burning interpulse stage, as 
implicitly assumed by the above core-mass - luminosity relationships) would not change the above conclusion 
that 56 Peg A cannot be a TP-AGB star, since these dips are about 0.3 dex fainter in logL (see Wagenhuber 
& Groenewegen 15 ). 

The conclusion that 56 Peg A cannot be a TP-AGB star plays a central role in the difficulty of finding for 
the 56 Peg system a suitable evolutionary channel which is consistent with the component masses derived 
from the spectroscopic orbit, as discussed in the next section. 

Difficulty with the masses under the corotation hypothesis 

The very small velocity semi-amplitude observed by Griffin (a mere 1.5 km s _1 ) yields a mass function 
of only 3.73 10 -5 M 0 . Under the assumption that the giant star’s rotation is synchronized with the orbital 
motion, and that both motions are in the same plane, Griffin uses the rotational velocity of 4 km s _1 (as 
derived by de Medeiros & Mayor 16 ) to obtain sin* ~ 0.22 (i.e., i ~ 13°). The Hipparcos astrometric data 
are unfortunately of no help to corroborate this value of the orbital inclination, as the sought astrometric 
orbit should have a± = 1 mas, just beyond the Hipparcos detection capabilities (especially when the orbital 
period is as short as 111 d; see Jancart et al. 17 ). The masses of the components are then constrained by the 
relation 

M 3 

---r ~ 0.0034 when i = 13°. (1) 

(Mi +M 2 ) 

From this relationship, Griffin derives a mass ratio q = Mi/M 2 of the order of 10 (where Mi and M 2 
are the masses of the giant component and of its companion, respectively), and lists pairs (Mi, M 2 ) = 
(2, 0.26), (4,0.40), and (6, 0.52) M 0 as examples of possible solutions. The author also notes that these 
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figures rely heavily on the above assumptions, but does not feel the need to reject them as they are indeed 
required for the system to remain detached in its present state. 

However, potential difficulties with the former evolution of the system inspired the authors of the present 
note to look more deeply into this result from a stellar evolution standpoint. More specifically, the possible 
masses derived from the mass function for the two components (under the assumption of i = 13°) raise 
several questions regarding the past history of the system, especially when considering the fact that the 
location of 56 Peg A in the H-R diagram forbids it from being a TP-AGB star, as discussed earlier. The 
very barium nature of 56 Peg A then implies that it has been polluted through mass transfer from a former 
TP-AGB companion. The present companion star has then to be a carbon-oxygen white dwarf (CO WD). 

This restricts the companion’s mass range to M 2 ~ 0.45 Mg. Therefore, among the three pairs of masses 
suggested by Griffin, only the high mass range is in fact relevant to 56 Peg. A semi-empirical initial final 
mass relation for WDs may then be used to further constrain M 2 . 

Taken together with the most recent initial-final mass relation of Weidemann 18 , expression m seems to 
require that the mass of the giant has been substantially altered by the mass transfer. Even so, however, 
there is hardly a way towards an acceptable solution. 

As an illustrative example, assume M 2 ,f = Myv d = 0.68 M 0 . Then from Weidemann relation, M- 2 a = 
3 Mq and from the mass function, Mi.f = 8.9 Mg. The former primary has thus lost in total 3 — 0.68 = 
2.32 Mq. Even assuming that all this mass was dumped on the companion yields an initial mass of Mpj = 
8.9 — 2.32 = 6.58 M 0 for the progenitor of the current giant star. But this mass is still larger than 
3 Mg, the estimated initial mass of the former primary. This solution is thus not acceptable, as it is 
obviously incompatible with the constraint that the WD progenitor must have been initially the more 
massive component. It is in fact possible to demonstrate graphically that there are no solutions for the 
whole range of relevant masses (the method is presented in what follows, where it is used with even more 
liberal assumptions on the theoretical final masses). Then we are left with the hypothesis that the initial- 
final mass relationship does not apply here. It is not surprising indeed, that in a mass-transfer system, a 
standard single-star M~Mf relation does not hold. We nevertheless note that the precursor of the current 
WD had to pass through the TP-AGB phase, in order to produce the heavy elements now observed in the 
atmosphere of the barium star. This again constrains the WD mass, as it is required that this mass be no 
smaller than the CO core mass at the first thermal pulse on the AGB, A-f c .iTp(M 2 .i). As indicated in the 
notation, this quantity is a function of the initial mass of the star. This relationship can be derived from 
stellar evolution theory and various sources give similar results - see e.g., Fig. 1 of Weidemann 18 . To fix the 
ideas, Fig. presents M Cj itp as a function of M 2 ,; from Blocker 19 (it is a compilation of his own data with 
those of Lattanzio 20 ) and from calculations by Herwig & Blocker (priv. comm, to Weidemann 18 ). Other 
theoretical results fall generally in between these two lines, so the range displayed on Fig. Ogives an estimate 
of the theoretical uncertainty. For our purpose the lowest theoretical predictions of M Cj itp for any given 
M 2 .i are relevant, as they are more liberal - in the sense that they allow a greater range of possible final 
solutions. If we conclude that it is hard to find a solution with Blocker’s curve, it will be even harder for 
any other existing relation. So, in what follows, Blocker’s relation 19 will be used. 

As can be seen, the absolute lower limit on the companion’s mass M 2 is now raised to ~ 0.52 M 0 and this 

value increases further with increasing M-iy This makes the present barium star quite massive (Mi ~ 6 M 0 
when the nominal inclination derived by Griffin is adopted). In fact, it becomes impossible to account for the 
formation of such a system in the framework of current stellar evolution theory. This difficulty is illustrated 
in Fig. 21 presenting mass-mass diagrams for the components of the 56 Peg system, for four evolutionary 
scenarios. In this figure, the horizontal axis corresponds to the mass Mi of the now-bariunr star, and the 
vertical axis to the mass M 2 of the now-WD. The actual evolution of the system would be represented by 
a curve joining points corresponding to initial and final configurations. As the exact course of the 56 Peg 
evolution is unknown, we consider two extreme scenarios, which should encompass the actual evolution: no 
mass transfer, i.e., only mass loss from M 2 (panels a and c) or conservative mass transfer from M 2 to Mi, 
i.e., no systemic mass loss (panels b and d). With the assumed mass-loss/mass-transfer modes, the system 
evolves along straight lines in the (Mi, M 2 ) diagram: vertical for pure mass loss, or inclined at —45° for 
conservative mass transfer. To guide the eye, the course of evolution is exemplified by several thin dotted 
lines in each panel. 

Adopting a given value for the initial mass ratio qi = M 1 J/M 24 , one may then represent the initial states 
of the system for any given qi as straight lines of appropriate inclination going through (0, 0). These are 
shown as thick lines, for qi = 1 in panels a and b, and for qi = 0.5 in panels c and d. Note that qi = 1 is the 
limiting case, as initially M 2 had to be at least slightly larger than Mi in order for the primary component 
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to evolve faster. 

Two independent constraints - which should be consistent with each other - are available to characterize 
the final (i.e., present) state of the binary system. First, the mass function expressed by Eq. 0 corresponds 
to the masses consistent with the observed radial-velocity curve and the adopted orbital inclination. This 
relationship between Mqf and M 2 y is shown in each panel as a thin solid line, accompanied by thin dotted 
lines representing the la errors (which combine the uncertainties on the mass function and on the inclination). 
Second, the theoretical relation ~M c ,itp (M 2) i) of Blocker (Fig. 0) provides another estimate of M 2j f (thick 
broken line) that connects to the initial value M 2i ; through the thin dotted lines describing the mass- 
loss/mass-transfer evolution. Only points at or above this line represent possible final configurations. Note 
that the transformation of this theoretical final locus due to the changes in assumptions is a simple stretching 
and squeezing along the horizontal axis. 

From Fig. 0 it is obvious that a marginal agreement between the observed and theoretical final masses 
is possible only in the most extreme case (qt = 1, conservative mass transfer: panel b of Fig. 0 ). But 
conservative mass transfer originating from the convective envelope of a TP-AGB star (the progenitor of the 
current WD) does not seem very plausible, as such a mass-transfer mode would rapidly suffer from dynamical 
instabilities. Another difficulty with this solution is the rather large mass implied for the current giant star 
(Mi > 5 M 0 ). According to the stellar-evolution grid of Lejeune & Schaerer 11 displayed in Fig. 0 such 
massive giant stars (of solar or subsolar metallicity) have absolute magnitudes well in excess of the value 
inferred for 56 Peg A. 

In fact, these tracks impose Mi < 4.5 M 0 for 56 Peg A. Such a limit on Mi can be inferred by inspection 
of the local My minima at the base of the giant branch: tracks for any star more massive than 4.5 M 0 will 
never reach as low a luminosity as My ~ —1.3 on the giant branch, irrespective of the metallicity. Thus 
the high mass solutions inferred from the spectroscopic orbit may be excluded right away from luminosity 
considerations. 

Considering the error bars on 56 Peg photometry and metallicity, the mass range 2-4 M 0 appears likely 
for the giant star (Fig. 0 . Then from Eq. 0 (corresponding to an orbital inclination of i = 13°), it follows 
that M 2 lies in the range 0.26-0.40 M 0 - much too small for a CO WD. 

Towards a new model 

How is it then possible to reconcile theory and observations? The V sin i estimate of de Medeiros & 
Mayor 16 seems beyond doubt, even though it involves some assumptions on the macroturbulent contribution 
to the observed line width. An independent determination by Gray 21 , using the Fourier transform technique 
which makes it possible to separate the macroturbulent from the rotational contribution, yields 3.9 km s _1 , 
very close to the 4.4 km s _1 value of de Medeiros & Mayor 16 . 

One simple possibility would then be to allow sin* outside the ltr range. In this case lower sin* values 
are needed: the more edge-on the system is, the higher the companion mass becomes and the higher the 
observational (Myf,M 2i f) locus lies on the mass-mass diagram. 

Dropping the assumption of corotation, involved in deriving the value of sin*, appears to offer very 
promising prospects, as we now show. This would improve the situation only if the giant rotates faster 
than corotation, implying a larger equatorial velocity. The inclination derived from the same value of the 
projected rotational velocity as adopted by Griffin would then become smaller. The observational final locus 
of 56 Peg on the (Mi, M 2 ) diagram would then come closer to our theoretical limit. To fix the ideas, for 
the qi = 0.5 and pure mass-loss case (case c), consistency between the mass-function and M c .iTp(M 2i i) loci 
can be achieved if * = 5°. This in turn implies an equatorial rotation velocity of about 48 km s -1 ! The 
only way to reach such a fast rotation is by spin accretion during the mass-transfer episode 22,23 . The spin 
accretion had to occur recently, because the dynamo associated with such a fast rotation and the resulting 
magnetic braking would slow the rotation down on a time scale of the order of 10 8 yr (see the paper by 
Theuns et al. 22 and references therein). One similar case is known among barium stars, namely HD 165141, 
a rapidly-rotating {V sin* = 14 km s _1 ) star, also classified 22,24,25 as RS CVn despite its long orbital period 
of 5200 d. The WD companion is the warmest known (T e ff ~ 35000 K) among barium-star companions, and 
its rather short cooling time scale (about 10' y) clearly hints at a recent mass (and spin) accretion event. 
One notes a large difference in orbital period between this system and 56 Peg, which arguably somewhat 
limits the resemblance. Spin accretion in HD 165141 is thought to occur through wind and accretion disc 
formation while the separation of components in 56 Peg is comparable to (or even smaller than) a typical 
AGB giant radius. Thus a mass transfer mode other than wind accretion had to be involved in the 56 Peg 
system. But it is worth noting that fast rotation (up to 100 km s -1 !) is in fact common among at least one 
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class of related systems, namely d’-type symbiotics 26 . These objects also do exhibit some degree of barium 
enrichment. The only d’-type symbiotic with known orbital period, V417 Cen, has V sini = 70 km s -1 and 
P = 247 d , which is much closer to the orbital period value of 56 Peg. That proves that spin accretion is 
somehow possible (and efficient!) even in systems with such relatively short orbital periods. 

So, could 56 Peg be a similar case? Very likely, as 

(i) the strong UV continuum 27,28 has been ascribed to a warm WD of T e s = 32000 ± 4000 K and M 2 = 
1 M 0 . Sweeney 29 predicts a cooling time of less than 50 10 * 6 y for a WD of that temperature. This 
cooling time is shorter than the magnetic braking time scale of the giant, so that the latter is still 
rotating fast. 

(ii) many different observers have reported signatures of chromospheric activity (see Griffin 1 and references 
therein). 

To mention just a few, Schindler et al. 27 , quoting Linsky et al. 30 , indicate that ‘56 Peg exhibits an unusual 
Ca II K line emission profile for its MK spectral type (...) The Ca II emission more nearly resembles that 
seen in solar plages.’ Eaton 31 has remarked that the Ha line profile is much shallower than the same line 
in inactive stars, thus suggesting that Ha has been filled in by emission. Similar indications of emission 
exist for the ultraviolet Mg II h & k doublet and the Na I D line (see Griffin 1 and references therein). The 
question nevertheless arises whether this activity is triggered by fast rotation or by mass transfer in the 
binary system. Several arguments support the hypothesis of a fast rotator. The first is the constancy of 
the Ca II H & K emission profiles over several years 27 , whereas a mass-transfer origin for these lines would 
rather cause some variability, associated with irregularities in the mass transfer.f Second, the bright UV 
emission lines of C II, C III, C IV, He II, Si II, Si III, Si IV, N III and N V observed by Schindler et al. 27 
along with the X-ray flux are typical signatures of a hot corona associated with magnetic heating and fast 
rotation. Schindler et al. 27 felt obliged to reject this obvious explanation, since the observed V sini is low. 
But the present paper precisely argues that these difficulties, as well as those raised by Griffin’s orbital 
elements, vanish if one assumes that the star is a fast rotator. Very interestingly, Schindler et al. 28 argue 
that 56 Peg deviates from the rotation-activity (Lx versus V sini) correlation of Pallavicini et al. 34 . Now 
it becomes obvious that it is the small sini value that drags this object that far to the left of the graph. 
Using corotation velocity (i.e. V = 18 km s _1 ) instead of V sin* improves the agreement, but 56 Peg is still 
an outlier (and would be even more so if one could resolve V from sin* for the other objects). However, 
adopting V ~ 50 km s _1 would restore 56 Peg right within the expected locus of RS CVn binaries, as far 
as its X-ray flux is concerned! This is shown in Fig. [3 which compares the X-ray flux of 56 Peg (corrected 
according to the new distance estimate from Hipparcos: 165 pc instead of 215 pc as estimated by Schindler 
et al. 27 ) to that of late-type (including active) stars. 

Given the small orbital separation, one may still wonder whether there should not be as well some 
contribution to the X-ray flux coming from accretion, as advocated by Schindler et al. 27 and Dominy & 
Lambert 35 . The major evidence suggesting that there may be a link between mass transfer and X-ray 
emission in active binaries is the correlation between the X-ray luminosity and the Roche-lobe filling fraction 
T 2 for RS CVn systems noted by Welty & Ramsey 36 . Singh et al. 37 re-examined this issue using samples of 
RS CVn and Algol binaries, and confirm the weak correlation found earlier by Welty & Ramsey 36 . However, 
Singh et al. “ regard this as a rather weak argument, because both the X-ray luminosities and the Roche lobe 
filling fractions are themselves correlated with the stellar radii in the sample of RS CVn binaries, and thus 
regard this correlation as merely a by-product of the inherent size dependence of these quantities. ” 

Actually, despite the short orbital period, the filling factor of 56 Peg A is not that close to unity,* and 
therefore mass transfer to the companion occurs only through the accretion of the giant wind, which is quite 
moderate according to Reimers’ mass-loss law: a mere 2.7 10~ 9 M 0 y _1 (adopting for the giant = 3 M 0 , 
R\ = 40 Rq and log L/Lq = 2.7). Schindler et al. 27 and Dominy & Lambert 35 showed that such a mass 
loss rate could indeed account for the observed X-ray flux in a 100 d system, using the Hoyle & Lyttleton 39 

t Surprisingly, the Mg II h & k doublet, which is supposed to form at somewhat higher chromospheric levels than the Call H 

& K, according to Mullan &: Stencel 32 , is strongly variable 27,32,33 . This variability cannot be ascribed to orbital modulation, 
since the system is seen nearly pole-on in the present model. Variations in the outflowing wind of the K giant, where the Mg II 
doublet supposedly forms, represent therefore an interesting possibility, as advocated by Mullan & Stencel 32 and Schindler et 
al. 27 . 

tAdopting i = 5°, Mi = 3 Mg, M 2 = 0.96 Mg, Griffin’s orbital elements yield an orbital separation A = 152 Rg and 
a Roche radius around the giant component of 73 Rg, well in excess of the 40 Rg representing the stellar radius itself, and 
translating into T 2 = 0.55. One may note, however, that rotation faster than orbital motion complicates somewhat the picture: 
the usual Roche geometry does not strictly apply, hence the effective T 2 will be a bit larger 38 . 
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accretion cross section. They failed, however, to recognize that this accretion cross section overestimates by 
about one order of magnitude the actual cross section, as inferred from hydrodynamical simulations 22,40,41 . 
This constitutes a further argument against a strong contribution from mass accretion to the X-ray flux. 

Symbiotic systems provide another interesting sample to which 56 Peg may be compared. One of the 
defining properties of symbiotic systems is to host a hot compact star, mostly a WD, accreting matter from 
a giant companion 42 . The WD is heated either directly by accretion, or indirectly by nuclear burning fueled 
by accretion (see Jorissen et al . 43 and references therein). Symbiotic stars are X-ray sources 44 , at the level 
10 3 ° to 10 33 erg s _1 , to be compared with 2.3 10 31 erg s _1 for 56 Peg 2 '’ 45,46 (corrected for the new Hipparcos 
distance). The physical process emitting these X-rays in symbiotic stars is still debated. Direct evidence for 
accretion disks, in the form of continuum flickering and far UV continuum (and, by extension, X-rays) is not 
usually found in symbiotic stars 47-49 , unlike the situation prevailing in cataclysmic variables and low-mass 
X-ray binaries 47,50 . Other mechanisms were therefore advocated to account for the X-ray emission from 
symbiotics 44 , like thermal radiation from the hot component in the case of supersoft X-ray sources, or the 
shock forming in the collision region between the winds from the hot and cool components. Interestingly, 
Soker 51 even suggests that the X-ray flux from symbiotic stars exclusively arises from the fast rotation of the 
cool component spun up by wind accretion from the former AGB component (now a WD)! If that hypothesis 
is correct, the properties of X-rays from symbiotics would be undistinguishable from those of active binaries. 

To conclude, it is thus unlikely that mass transfer possibly occurring in the 56 Peg system contributes 
significantly to the observed X-ray flux. The X-rays from 56 Peg share many properties of active binaries. 
On top of the arguments mentioned above, the ROSAT HR1 hardness ratio for 56 Peg, as given by Voges 
et al . 45 (HR1 = 0.68, where HR1 — (H — S)/(H + S ) and H and S are the ROSAT PSPC count rates in 
the 0.5 - 2.0 keV and 0.1 0.4 keV ranges, respectively) does lie in the range seen in RS CVn stars taken 

from the list of Singh et al. 37 . 

Another argument supporting the presence of rotation-driven activity in 56 Peg is provided by the surface 
flux Sciv °f the C IV A155 nm line, which falls exactly on the relation S'civ — P I0 t derived by Gunn et al . 52 
for active binary stars. Adopting for 56 Peg A a radius of 40 R 0 and a distance of 165 pc, the observed 27 
C IV flux [1.2 10 -12 erg cm -2 s _1 ] translates into a surface flux of 4.0 10 4 erg cm -2 s -1 , as compared to 
5.1 10 4 erg cm -2 s _1 from the period - surface-flux relationship of Gunn et al . 52 for an estimated rotation 
period of P ro t = 42 d, adopting Vot = 48 km s _1 . Similarly, the Sx/Sboi and Sciv/Sboi ratios for 56 Peg 
fall on the tight relation obtained by Dempsey et al . 53 for RS CVn systems. 

Conclusions 

The following coherent picture thus emerges for the 56 Peg system: 

• Evolutionary tracks combined with Tycho-2 photometry and Hipparcos parallax (< V >= 4.76, < 
B — V >= 1.28, w = 6.07 ± 0.67 mas, My = —1.32)J)g' 2 g) imply that the giant component in the 
56 Peg system has a mass of the order of 2 to 4 M@, for metallicities solar or slightly subsolar (—0.4 < 
[Fe/H] < 0.). At solar metallicity, the Lejeune & Schaerer 11 tracks predict M = 3 M 0 for our adopted 
56 Peg A values of T e g = 4200 K and logL = 2.7, corresponding to Mboi = —2.0. 

• Whatever its metallicity in the above range, 56 Peg A is not luminous enough to be on the TP-AGB. 
The barium overabundance observed in the atmosphere of 56 Peg A must therefore be ascribed to mass 
transfer in the binary system. 

• The WD companion of temperature 32000 K (inferred from a fit to the observed UV continuum) 
implies that the WD is young (with a cooling time scale of the order of 50 10 6 y). 56 Peg A might 
therefore have been spun up during the mass transfer process (as d’-type symbiotics of similar orbital 
periods, thus representing yet another example of the class of WIRRING systems defined by Jeffries & 
Stevens 23 ) and not yet slowed down by magnetic braking, which occurs on a time scale longer than the 
WD cooling time scale. In this scenario, high rotational velocity is achieved after and independently 
of tidal orbital circularization. 

• 56 Peg A is therefore an active star, as evidenced by many emission lines in the optical (Ha, Na D, 
Call H & K) and ultraviolet. Diagnostics like Sx/Sboi, Scrv/Sboi (where S denotes the surface flux) 
and Lx versus Vot all fall along the relations expected for active stars of the RS CVn type. The X-ray 
flux is therefore attributed to the activity of the giant star rather than to mass transfer in the binary 
system. 


6 


• The strong variability of the Mg II h & k lines has been ascribed to variability in the wind from the K 
giant. The above stellar parameters predict 3 10 ~ 9 M 0 y -1 from the Reimers mass-loss law. 

• All activity indicators are consistent with a rotational velocity of 30-50 km s , implying an orbital 
inclination i ~ 5-8°. 

• The orbital mass function combined with the giant mass range (2-4 M 0 ) then implies M 2 = 0.75 to 
1.15 M 0 , now consistent with it being a CO WD. 

• If the system is seen nearly pole-on, Griffin’s orbital elements then yield an orbital separation A = 
152 R 0 (for Mi = 3 M 0 and M 2 = 0.96 M 0 ), implying a Roche radius around the giant component 
of 73 R 0 , well in excess of the 40 R 0 representing the stellar radius itself. The short orbital period 
(111 d) is not at all incompatible with the system being currently detached. 

• The assumed inclination, * = 5°, has an a priori probability of one part in 260. Admittedly, it is 
a small value, but it does not make the proposed solution impossible. In the light of the presented 
evidence, and bearing in mind the opening quotation from Sir Arthur Conan Doyle, the authors feel 
encouraged enough to submit their reasoning under the reader’s judgement. 
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P (d) 


Figure 1: The eccentricity - period diagram of barium stars, with 56 Peg represented as an open square. 
Data for the other stars are from Jorissen et al. 54 . 
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B-V B-V 


Figure 2: Position of 56 Peg A from Hipparcos and Tycho-2 photometry (black dot within an error box; 
note that the formal error in My is completely dominated by the parallax uncertainity) in the observational 
H-R diagram, plotted against Geneva theoretical tracks 11 , of various metallicities as indicated in the figure. 
Labels indicate initial masses of the stars. The thick oblique line protruding from the observed location of 
56 Peg A corresponds to the de-reddening vector for Ay = 0.15 mag, as derived from the extinction model 
of Arenou et al. 55 . 
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Figure 3: Mass of the CO core at the first thermal pulse, M Ci itp, as a function of initial mass M 2.0 according 
to Blocker 19 - line with filled squares - and to Herwig & Blocker (priv. comm, to Weidemann 18 ) - line with 
rhombs. 
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Mj Mj 

Figure 4: Mass-mass diagrams describing the evolution of the 56 Peg system. The horizontal and vertical 
axes correspond to the masses of the current barium star (Mi) and of the current WD (M 2 ), respectively. 
Each panel corresponds to a specific evolutionary scenario: (a) qi = = 1, mass loss only (no mass 

exchange); (b) qt = 1, conservative mass transfer; (c) qi = 0.5, mass loss only; (d) qi = 0.5, conservative 
mass transfer. See text for details. 
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Figure 5: X-ray luminosities vs. projected rotational velocities for 56 Peg (encircled crosses, for V sin* and 
for two values of V : assuming orbital corotation or fast rotation) and stars of various spectral types and 
luminosity classes detected by the Einstein Observatory. Different symbols indicate different luminosity 
classes. Adapted from Pallavicini et al. 34 . 
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